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Abstract 

Strain rate and temperature dependence of the tensile strength of single-wall 
carbon nanotubes has been investigated with molecular dynamics simulations. 
The tensile failure or yield strain is found to be strongly dependent on the tem- 
perature and strain rate. A transition state theory based predictive model is 
developed for the tensile failure of nanotubes. Based on the parameters fitted 
from high-strain rate and temperature dependent molecular dynamics simu- 
lations, the model predicts that a defect free //m long single-wall nanotube at 
300K, stretched with a strain rate of 1%/hour, fails at about 9 ± 1% tensile 
strain. This is in good agreement with recent experimental findings. 
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Carbon nanotubes (CNTs) have very high strength and stiffness under axial stresses, 
but are very flexible under non-axial stresses. The theoretical and experimental values of 
the Young's modulus are found to be around 1 TPa. For large tensile stresses, the atomistic 
molecular dynamics (MD) simulations performed at very high strain rates (limited by time- 
scale of the phenomenon that can be simulated in MD) show tensile failure strains to be as 
high as 3O%.0 Whereas, recent experiments report much smaller tensile failure strains in a 
variety of scenarios. Walters et al@ have reported a failure strain of 5.8% on single walled 
CNT (SWCNT) ropes. Yu et all have reported a similar value of maximum 5.3% strain 
for failure at room temperature, and have measured strain for failure of multi-walled CNTs 
(MWCNT)i to be as high as 12%. 

Earlier theoretical and numerical simulation studies have shown that, under large tensile 
strains, Stone-Wales(SW) bond rotations result in the formation of pentagon-heptagon pair 
(5775) defects on the nanotubes (Figure 1). The formation energy of such defects decreases 
with the increasing tensile strains and these defects are energetically favorable at tensile 
strains larger than 5%,i However, the large energy barriers (of about 8-9 eV) to the forma- 
tion of SW bond rotations in static calculations^ do not readily explain the experimentally 
measured low tensile strength of SWCNT ropes or MWCNTs. The high strain rates used in 
the MD simulations, on the other hand, are generally of the order of pico second' 1 or faster, 
which are unrealistic, as compared to the experimental strain rates of minute' 1 or slower. 
A clear deficiency thus exists in a direct comparison of MD simulated tensile strength with 
what is measured in experiments. It is not clear if the differences are due to the real physical 
and chemical mechanistic reasons, or due to the limitations of the extremely fast strain rates 
used in the atomistic MD simulations. 

A typical time step in MD simulations can not exceed a fraction of the vibrational 
periods of atoms or molecules under investigations, and is generally of the order of few femto- 
seconds or less. Hundreds of millions of time steps, therefore, are needed to simulate typical 
nanosecond dynamic processes in atomic systems. For thousands of atoms, interacting with 
complex many-body force field functions, the MD simulation is computationally intensive 
and is not always accessible for routine investigations. One way to accelerate the kinetic 
processes, that occur during MD simulations, is to increase the temperature. The transition 
time, for an activated process, to go from one state to another is given by the Arrhenious 
relation t = l e E ^/ k BT ^ w h ere j s the activation energy and v is the effective vibration or 
attempt frequency of the transition. At high temperatures the larger kinetic energy increases 
the rate of the process to overcome the barriers and the transition time is thus shortened. If 
the tensile yielding of a material can be described by a series of activated processes of defect 
formations. It may be possible to describe the tensile yielding or failure of CNTs with an 
Arrhenious type model incorporating both short time and high temperature defect formation 
behavior, as well as long time and low temperature behavior. We have used temperature and 
strain rate dependent (with in MD time-scale) MD simulations to propose and parameterize 
a predictive model for tensile yield strain of micron long SWCNTs at experimentally feasible 
strain rates and temperatures. 

The MD simulations involve continuous straining of nanotubes by moving the end atoms 
under given strain conditions, such as under tensile and compression strain, and letting the 
nanotube relax during the dynamics at constant temperatures and strain rates. Tersoff- 
Brenner potentialiJi is used to model Carbon-Carbon (C-C) interactions. A 60A long (10,0) 



2 



nanotube with 600 atoms is used for the simulation studies, and the CNT is strained at 
rates varying from 10 _6 |?s _1 to lO^ps^ 1 and at temperatures ranging from 300K to 2400K. 

In Figure 2, the strain energy per atom for different temperatures and strain rates, as 
a function of tensile strain, is plotted. The abrupt deviation from the elastic behavior, at 
high strain values (and we define such strains as yield strains or failure strains), is due to 
the formation of SW 5775 defects followed by necking, thinning and eventual breaking of 
the nanotubes. For strain values higher than the yield strain, the perfect hexagonal config- 
urations of the carbon atoms in the CNTs are destroyed. This is in qualitative agreement 
with the previous MD simulations. 00 Figure 2 also shows that the yield strain (and the 
corresponding tensile strength) decreases at higher temperatures and at slower strain rates. 
Initial MD studies by Yakobson et. al.0 showed similar behaviors but with very high yield 
strains due to the very high strain rates used in their simulations. Our simulations with a 
larger variation of strain rates and temperatures show the changes in the yield strain (or 
tensile strength) in a wider range. For example, at slow strain rate (3 x 10~ 5 /ps) and high 
temperature of 2400K, the nanotubes yield at about 6% in our simulations. The yield strain 
as a function of strain rate, at different temperatures, is shown in Fig. 3. The linear depen- 
dence at each temperature indicates that the Arrhenious behavior is valid and the tensile 
yielding of nanotubes can be described by a theory of activated processes, i.e., the transi- 
tion of the system from the initial pre-yielding state to the final post-yielding state occurs 
through a series of activated processes with an effective barrier defining the characteristics 
of the whole phenomenon. 

According to the Eyring's theory the Arrhenius formula for the strained system can be 
modified asQ 

1 1 E V -VK* 

t = e , (1) 

where e is local strain, n site is number of sites available for state transition, K is force 
constant and V is the activation volume. Thus the strain rate required for a transition at a 
certain strain is: 

1 11 E v -VK t 

— = — e k B T , (2) 

where e is strain rate and eo is a constant related with vibration or attempt frequency. 
The yield strain as a function of temperature and strain rate, thus, can be expressed by 
inverting Eq. (2) as, 

eY -VK + VK ln ^J- (3) 

For nanotubes under tensile strain, we have observed that the yielding process induces 
a sequence of S-W type bond rotations within the connected local regions. These multiple 
transitions lead to a collective kinetic activation mechanism of the failure of CNTs strained 
beyond elastic limit. The combined rate for N multiple transitions thus can be of the form 
4t, where | is the rate for a single transition and for simplicity we have assumed that the 
rate for each successive transition can be expressed by an averaged value of 1/t. An averaged 
effective activation energy E v for the multiprocess thus replaces E v for the single process, 
and Eq. (3), can be modified as: 
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To fit the data from the simulations, Eq. (4) can be rewritten as: 

ey = hi — H In 

K VK VK y NAe step " VK Ae step 



= A{T) + B{T)\n-—, (5) 

^step 

Mv _ k BT ( n sit< ,e Q \ 1 d/t->\ _ k B T 



where A(T) and B(T) are defined as A(T) = f% - w ln ( jVAel" ) and B ( T ) 
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Ae s t ep is the change of strain at each step used in simulations. In Fig. 4, the functional 
dependence of A(T) and B(T) is obtained over the entire temperature range by fitting the 
linear dependence to the simulation data. The linear dependence of B{T) on the temperature 
indicates that the transition state theory (within Eyring's model) is valid for the yielding 
of the strained CNTs. From the dependence of B(T) on T, we get VK = 18.04 eV. For 
one-dimensional CNTs, with Young's modulus giving the force constant to be 1 TPa, the 
activation volume comes out to be 2.88A . This corresponds to a typical atomic volume of 
a carbon atom within a nanotube. The fitting of the coefficients in A(T) gives = 0.20 
and ln( ^"£7 ) = 3.66 x 10~ 5 . These values correspond to the average kinetic activation 
energy E v of 3.6 eV for the successive multiple S-W bond rotations. In our simulations, 
with n site ~ 600 and Ae step = 0.0025 or 0.25%, ^ comes out to be about 8 xl0~ 3 ps^ 1 . ey 
is not dependent on Ae s t ep , as the effects cancel out in A(T) and B(T). There is always a 



possibility for the inverse transitions which have rates scaled as e k B T e, and are omitted as 



>> 1 in our study. 
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The yield strain of the CNTs, under tensile stress with experimentally feasible strain 
rates, can be estimated or predicted from the above model. For example, at T = 300K and 
at a strain rate of about 1% per hour the yield or failure strain for a 6nm long (10,0) CNT 
comes out to be about 11 ± 1% as according to ey (%) = 20.21 + (0.21 ±0.02) lne {ps~ l ) from 
the data shown shown in Fig. 3. Furthermore, the effect of experimentally realistic length 
of CNTs on the predicted values of ey can be estimated from Eq. (4) as according to, 

Ae = "V^ ln ( n ^/ n °site), (6) 

where n s u e is linearly dependent on the length of CNTs. For a fim long (10,0) CNT, 
under the same conditions, the yield or failure strain lowers to about 9±1%. Of course, much 
more total strain energy is needed to yield a longer CNT than the energy needed to yield 
a shorter CNT. A larger diameter nanotube of the same length may have more activation 
sites n S ite for nucleating the defects leading to the yielding of the tube. However, N, the 
number of single processes involved in the yielding, may also increase and offset the effect. 
This will not be the case if individual single-wall nanotube diameter increases by orders of 
magnitude, which is highly unlikely. 

In Eq.(4), we note that the yield strain of CNTs is linearly dependent on the average 
kinetic activation energy. The static in-plane activation energy of the first S-W bond ro- 
tation, for several armchair and zigzag CNTs, were calculated at zero tension, by rotating 
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the C-C bond in the CNT surface plane and relaxing the system with zero temperature 
MD simulations. The results are shown in Figure 5. E u is found to be around 9.5 eV for 
large diameter nanotubes. This is about 1 eV different than the activation energy of 8.5 eV 
computed by Zhang et. al., for a (6,6) nanotube using static tight-binding simulation.! This 
means that the error in the computed averaged activation energy of the multiple kinetic pro- 
cesses leading to the yielding of the nanotube, due to the use of Tersoff-Brenner interactions 
(as opposed to a different and perhaps more accurate tight-binding or ab-initio interactions) 
could be estimated to be as high as 1 eV, and can give an estimation of an additional error- 
bar on our predicted value of the yield strain. According to Eq. (4) 1.0 eV error in the 
averaged E u corresponds to about 5% error in the yield strain. The predicted yield strain 
from the MD simulations using Tersoff-Brenner interactions can further be improved with 
more accurate tight-binding molecular dynamics simulations in future. Moreover, we note 
that the static activation energy of the first SW bond rotation is more than twice higher 
than the averaged kinetic activation energy of successive kinetic S-W bond rotations. This 
is natural, because the activation energy of the S-W bond rotation, in a strained tube, and 
in the presence of already existing S-W bond rotations or other defects, is expected to be 
much lower than the activation energy of the first S-W bond rotation or other defects in an 
ideal perfect tube. 

Additionally, the plot in Fig. 5 also shows that a small diameter zigzag nanotubes has a 
smaller activation energy than that of a large diameter zigzag nanotube. No such difference 
is obvious for the armchair nanotubes. It is expected that the yielding behavior will show a 
minor diameter dependence for the zigzag nanotubes and no such dependence in the armchair 
nanotubes. This is shown in the inset in Fig. 5. This also shows that the averaged kinetic 
activation energy for successive multiple kinetic S-W defect formation processes is a very 
important parameter in our model and could be improved upon in higher fidelity dynamic 
or static simulations in future. 

Above analysis shows that the yield strain of a tensile strained CNT is strongly dependent 
on the temperature and the applied strain rate. Figures 3 and Eq.(4) show that slower 
strain rates at lower temperatures can be equivalent to faster strain rates but also at higher 
temperatures. Both the conditions increase the probability of overcoming the activation 
barriers during the dynamic processes. From Eq.(4), this equivalence can be defined as 

( Jl^n = (WL)* (7) 

and the processes can thus be accelerated by elevating the simulation temperatures. In 
a different context, Voter et. al. haveused the high temperature accelerated dynamics 
to study the diffusion on solid surf aces, EJ'Ej which shows that more realistic experimental 
deposition rates can be reproduced through this scheme. Of course the dynamics on CNTs 
is much simpler compared with other materials, as CNTs are quasi-one dimensional and 
have simple structures. 

For comparison, the strain rate and temperature dependence of the yielding of CNTs 
under the compressive stress is also studied using the MD simulations with Tersoff-Brenner 
potentials. Sideways buckling of the compressed nanotube at lower temperatures, similar 
to Yakobson et al's MD study at zero temperature,0 and a collapsing behavior at higher 
temperature driven by diamond-like bond formations, as noted by Srivastava et. al. in 
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quantum molecular dynamics simulations were observed.li.il Additionally, S-W type bond 
rotation defects were also observed in simulations at higher temperatures.^ These defects 
are found to be irreversible with the removal of the compression induced strain in the 
system. Our MD simulations show that the diamond-like bond formations and S-W bond 
rotation defects are strongly correlated with the local buckling instability of the compressed 
nanotubes, which tend to additionally reduce the barrier heights in the compression case. 
The complex coupling of the buckling formation with the activated processes, as described 
in this work, is currently under investigations and will be discussed separately. 

Finally, experiments show the tensile strength of 40-50 GPa for SWCNT ropes and 
MWCNTs, even though their yield strains are different: 5 — 6% for SWCNTs and 12% for 
MWCNTs. Eq. (4) shows that the yield strain ey is also a function of activation volume 
V and Young's modulus K. Since K is shown to be about 1 TPa for diverse SWCNTs 
and MWCNTs, the difference would be related to the activation volume. It is reasonable 
to expect a smaller activation volume on the outer shell of a MWCNT due to the presence 
of inner shell CNTs, and the reduced V will increase the value of ey for MWCNTs. These 
issues are currently under investigations and will be published in future. 
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FIG. 2. The strain energy per atom as function of tensile strain for CNT(10,0) at different 
temperatures and strain rates. 
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FIG. 3. Yield strain of (10,0) CNT as a function of strain rate at different temperatures under 
tensile stress. 
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FIG. 4. The fitting of A(T) and B(T) with data obtained from Figure 3. 
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FIG. 5. The static in-plane activation energy of Stone- Wales bond rotations as a function of 
diameters for zigzag and armchair CNT, using Tersoff-Brenner potentials. Inset (Figure 6 and 
Figure 7): The yield strain of (10,0), (20,0), (5,5) and (10,10) CNTs as a function of strain rate at 
T=2400K, from MD simulations. 
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FIG. 6. Inset of Figure 5: The yield strain of (10,0), (20,0) CNTs as a function of strain rate 
at T=2400K, from MD simulations. 
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FIG. 7. Inset of Figure 5: The yield strain of (5,5) and (10,10) CNTs as a function of strain 
rate at T=2400K, from MD simulations. 
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